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An upper limit to non-Newtonian attractive forces is obtained from the measurement of quantum 
states of neutrons in the Earth's gravitational field. This limit improves the existing contraints in 
the nanometer range. 
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I. INTRODUCTION 

According to the predictions of unified gauge theories, 
supersymmetry, supergravity, and string theory, there 
would exist a number of light and massless particles . 
An exchange of such particles between two bodies gives 
rise to an additional force. Additional fundamental forces 
at short distances were intensively studied, in particu- 
lar during last few years following the hypothesis about 
"large" supplementary spatial dimensions proposed by 
Arkami-Hamed, Dimopoulos, and Dvali building 
on earlier ideas in Refs 0, 0, 0, 0- For a review of 
theoretical works and recent experimental results, see 
IE H> ^> ^ ■ This hypothesis could be verified using 
neutrons because the absence of an electric charge allows 
one to strongly suppress the false electromagnetic effects 
[T^. It was noticed in [l^ that the measurement of the 
neutron quantum states in the Earth's gravitational field 
|l4j is sensitive to such extra forces in the sub-micrometer 
range. In the case of n = 3 extra dimensions, the char- 
acteristic range is just in the nanometer domain 0, 
which is accesible in this experiment. The first attempt 
to establish a model-dependent boundary in the range 
1-10 /im, was presented in ref. 16] . 

An effective gravitational interaction in presence of an 
additional Yukawa- type force is parametrized as: 



V^„{r) = G 



miTO2 



1 -I- ace" 



-r/X 



(1) 



Here, G is the Newtonian gravitational constant, mi and 
TO2 are interacting masses, r their relative distance, ac 
and A are strength and characteristic range of this hypo- 
thetical interaction. 

The experiment 14] consists in the measurement of the 
neutron flux through a slit between a horizontal mirror 
on bottom and a scatterer/absorber on top as a func- 
tion of the slit size A/i.This dependence is sensitive to 
the presence of quantum states of neutrons in the poten- 
tial well formed by the Earth's gravitational field and the 
mirror. In particular, the neutron flux was measured to 
be equal to zero within the experimental accuracy if the 



slit size A/i was smaller than the characteristic spatial 
size (a quasiclasssical turning point height) of the lowest 
quantum state of ~ 15/xm in this potential well. The 
neutron flux at the slit size A/i < 10/im was lower by 
at least a factor of 200 than that for the lowest quan- 
tum state(Aft, w 20/xm). If an additional short-range 
force of sufficiently high strength would act between neu- 
trons and the mirror then it would modify the quantum 
states parameters: an attractive force would "compress" 
the wave functions towards the mirror, while a repulsive 
force would shift them up. In this experiment, no devi- 
ation from the expected values was observed within the 
experimental accuracy. This accuracy is defined by the 
uncertainty in the slit size which can be conservatively 
estimated as « 30 % for the lowest quantum state ^j. 

The motion of neutrons in this system over the vertical 
axis z could be considered, within first and quite good 
approximation, as an one dimentional problem for which 
the mirror provides an infinitely high potential. The in- 
teraction between neutrons and the Earth is described by 
the first term in equation and can be approximated 
by the usual linear potential [r — R -\- z): 



V{z) ~ mgz 



(2) 



with g = GM/R^, R being the Earth's radius, M its 
mass, m the neutron mass. 

The second term in equation ^ introduces an addi- 
tional interaction. Due to the short range of this interac- 
tion, its main contribution is provided by the interaction 
of neutrons with a thin surface layer of the mirror and 
the scatterer. 

In this article, we estimate an upper limit on an addi- 
tional attractive short-range force, which could be estab- 
lished from this experiment in a model-independent way. 
We show that it could not be significantly improved in 
any more sophisticated model. 
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II. ATTRACTIVE INTERACTION 

Let us estimate firstly the interaction of neutrons with 
the mirror. If the mirror's density is constant and equals 
to then an additional potential of the interaction 
between neutrons and the mirror is given by: 



V'{z') = —Gaampn 



exp 



(3) 



The volume integral is calculated over the mirror bulk: 
—CO < x^y < GO, 2; < (in fact, over the neutron's 
vicinity with the size of the order of a few A due to the 
exponential convergency of these integrals). It can be 
calculated analytically for small A: 



V'{z) 



-Uoe 



-z/A 



(4) 



with Uq — 2'KGaGrapm)^ ■ 

The simplest upper limit on the strength of an ad- 
ditional interaction follows from the condition that this 
additional interaction does not create itself any bound 
state. It is known 15] that for an exponential attractive 
{aa > 0) potential Q this means that: 



UomX^ 



< 0.72. 



(5) 



This condition gives a boundary for an additional poten- 
tial strength: 



aa = 0.72- 



P 



h R 



TI" Pr. 



mgX^ mX A 



(6) 



p being the Earth's averaged density. In this experiment, 
both densities are close to each other p « p,n, therefore 
their ratio p/pm is close to 1. However an adequate choice 
of the mirror material (coating) would easily allow one 
to gain a factor of 3-5 in the sensitivity in future exper- 
iments. One obtains the following numerical boundary: 



15 I 



(7) 



Here, 1 pm is chosen as a natural scale for this experi- 
ment. This limit is presented in Figure ^ in comparison 
with the limits coming from the Casimir-like and van der 
Waals force measurement experiments One can note 
that, in the realistic case, one has to establish a condi- 
tion of non existence of an additional bound state for the 
sum of (0) and l@J but not for the interaction Q alone. 
A presence of the linear potential modifies slightly the 
critical value in jSJ. For instance, for A — 1pm it equals 
approximately to 1.0 and for A = O.l^m it equals to 0.74. 
For smaller A, this value is evident to tend to 0.72. It is 
possible to explain qualitatively why the strength of an 
additional interaction should be higher in presence of the 
TOg2;-potential than without it. When a bound state just 



appears, then its wave function is extremely spread. If 
one adds a supplementary "external" confining potential, 
it does not allow the wave function to be spreaded and 
thus one needs a stronger potential to create a bound 
state. The range of presented A is 1 nm - 10 pm. A 
deviation from a straight line in the solid curve at 1 nm 
is due to the finite range of increase of the mirror effec- 
tive nuclear potential (impurities on the surface and its 
roughness). The same effect at A « 10/im is due to an 
"interference" of the potentials ^ and ^ . 
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FIG. 1: The constraints on ac following from the experi- 
ment [T3 | (the solid line) in comparison with that from the 
measurement of the Casimir and the van der Waals forces 
(the dashed line) . The dotted line shows a limit which can be 
easily obtained by an improvement of this experiment. The 
solid horizontal line rep resents the limit established from the 
atomic experiment 

It is interesting to compare this analytical limit 
to an analogous expression obtained in ITol | for the 
Casimir force-like experiments. The simplest boundary 
a = a{X) following from these experiments is given by a 
formula: 



ac = Cc 



exp(rfo/A) 
A3 



(8) 



Here c?o is a- gap separation and Cc is a constant depend- 
ing on the geometry of different experiments. This func- 
tion increases exponentially when A tends to zero. This 
behaviour is clearly seen in Figure ^ In the experiment 
El we obtain 



A^ 



(9) 



and a„ increases only as 1/A''! This difference between 
and (j3 means that, in principle, for any Cc and C„, 
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one could find a domain of sufficiently small but exper- 
imentally observable A in which the limit obtained from 
our experiment is stronger than that obtained from a 
Casimir force-like experiment. 

Up to now, we have not discussed the scatterer. As 
we can see from Figure ^ a competitive limit could be 
established only for A much smaller than the scatterer 
roughness amplitude of the order of 1 /zm. Therefore an 
influence of the scatterer is negligible. 



III. REPULSIVE INTERACTION 

Unfortunately, this experiment docs not allow us to 
establish a competitive limit for a repulsive interaction. 
In this case, there could be no "additional" bound state. 
Here, instead of the condition of "non-existence" of a 
bound state, one could consider the critical slit size for 
which the first bound state appears in this system. Such 
an approach would be model-dependent due to uncer- 
tainties in the description of the interaction of neutrons 
with the scatterer. Nevertheless, it is possible to obtain 
a simple analytical expression for small A and to show 
explicitly a difference in sensitivity of this experiment to 
an attractive and to a repulsive additional interactions. 

Let us remind the reader that the quasiclassical ap- 
proximation for the potential composed by the linear 
term mgz and by the infinite mirror potential describes 
the problem very well. For instance, it gives the ground 
state energy, for which the quasiclassical approximation 
is usually expected to be not valid, within 1 % [TEL l2?il|. 
The Bohr-Sommerfeld quantification condition in this 
case assumes that the quasiclassical wave function 



(10) 



with p{z) — ^2m(E — mgz) {0 < z < b, b being the 
turning point) equals to zero at z = (at the inpenetrable 
mirror). This condition, for ri-th quantum state, means 

p{z)dz + - = T:{n + I), n = 0, 1, . . . (11) 



and the corresponding energies are: 

En = i^^y\mg^hY/Hn + 3/'i)y^ (12) 

With an additional short-range repulsive interaction 
V'{z) — J7o exp{— z/A}, the problem could be solved 
in an analogous way. One can choose such zq that 
A <C ^0 ^- For < 2; < zq, the mgz potential is negligi- 
ble and the problem can be solved analytically ll5j. For 
zq < z < b, one can use the quasiclassical approxima- 
tion. At z = zq the quasiclassical wave function should 
be equal to the exact one: 



*cx(2o) = C2 sin(fczo + So), 



(13) 



with the phase shift 
60 -4fcA 



C 



2 loiO 



(14) 



Here k = ^2mE/h'^, ^ = ^^mUoX^/K^ , C = 0.577 
is the Euler constant, /o(C) and ifo(C) are the modified 
Bessel functions of pure imaginary argument. This phase 
shift modifies the Bohr-Sommerfeld quantification condi- 
tion and introduces a shift 5En in the level's energy. If 
this shift is small with respect to the nonperturbated 
value, one obtains the following relation between the en- 
ergy shift 5 En and the phase shift 5q: 



2El 



mg'^ 



1/2 



= Sa. 



(15) 



For very strong repulsive potential, ^ 1 , the phase 
shift is equal to 



So 



-4/cAln 



(16) 



An experimental upper limit on the energy shift SEn 
would impose an upper limit on ac for a repulsive in- 
teraction: 



with 



or 



2UomX'^ 



Ao 



< exp(Ao/A). 

^ 5_En 
mg ' 



o-G < - 



h R 



TT mgX^ mX X 



exp(Ao/A). 



(17) 



(18) 



(19) 



Direct comparison of relation (|19|l to |5j and © shows 
that the limit (|19|) at small A is sufficiently less restrictive 
than that for an attractive one © due to the exponential 
factor. 

Note that the presented approach could be ap- 
plied as well for an attractive interaction V'{z) — 
— [/q exp{— z/A}. In this case, the phase shift is equal 
to 



-4fcA 



C 



2Jo(e) 



(20) 



Here Jo{£,) and iVo(^) are the Bessel functions. This ex- 
pression is maximal in the vicinity of the first Bessel func- 
tion's zero ^0 ~ 2.40 where it can be developped as 



1 



2 M^o) e-Co' 



(21) 



One thus obtains the relation between the strength of an 
attractive potential and the energy shift SEn < 0: 



Ji(^o) Ao' 



(22) 
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with Ao = Noi^o) « 0,51, Ji(Co) « 0,52. If one 

omits the second term in the latter equation, one obtains 
exactly the relations © and © representing a condition 
of non-existance of an additional bound state. 

Let us mention that this formula allows us to estimate 
the capability of this experiment to further improve an 
upper limit on aa in the case of an attractive interac- 
tion: in order to improve significantly (by an order of 
magnitude) the limit © in a given range of A it is nec- 
essary to determine experimentally and theoretically the 
parameters of quantum states with a precision Aq/A « 1. 



IV. OCCUPATION NUMBERS 

The considerations presented above are valid only if 
the neutron population in the lowest quantum state in 
such a system (with an additional interaction included) 
if sufficiently high to provide a measurable signal/noise 
ratio. As we have mentioned in the Introduction, the 
experiment [3 | would allow one to identify an additional 
quantum state if its occupation number would not be 
suppressed by more than a factor of 200 compared to 
that for other states. 

In order to calculate the occupation numbers, let us 
start with a general expression for the probability of a 
rapid transition from a state k with the wave function 
^'fe(a;) to a state n with the wave function $„(a;) which 
is given by a formula [2ll |: 



Wk- 



(23) 



For a few initial quantum states, the probability Wn is a 
sum (an integral for continius spectrum) over them: 



(24) 



with the occupation numbers fk of initial states. 

To obtain an analytical expression for the occupation 
numbers, let us consider a simplified model of a harmonic 
oscillator in a final state and a plane wave in an initial 
one. An explicite analytical shape of the final state wave 
function does not play a role (the only important parame- 
ter is its spatial size) and would not modify considerably 
the occupation numbers. The used wave functions are 
equal to 



5'n(a;) 



1 



1 



: exp 



2xq 



a;o 



and 



*fe(a;) 



(25) 



(26) 



Here Hn{z) are the Hermite polynomials and xq repre- 
sents a characteristic geometric size of the problem. 



If initial states are populated according to the Gaus- 
sian law with a chercteristic momentum fcp then 



/fc = exp ( --2 



(27) 



and all integrals 1)24(1 can be calculated analytically [2^ . 
The occupation numbers of the final states appear to be 
equal to: 



kpxo / (fcpxo)^ - 1 
v/1 + (koxo)^ \ikoXo)^ + l 



i/2 



XP, 



^{koXo)^ - 1 



(28) 



Here Pn{z) is the Legendre polinomial. For instance, for 
the lowest states with n = and n = 1: 



Wo = 



koxo 



y^l + (koXf))' 



(29) 



If koXo ^ 1 then the occupation numbers are approxi- 
mately equal for all states: 



1. 



(30) 



Let us apply these expressions to our problem. For 
the gravitational quantum states, Xq ~ 6/im; the verti- 
cal velocity distribution has a characteristic velocity of 
vq « 50 cm/s. For these states, fco^o ~ 50 ^ 1 and all 
states should have approximately the same occupation 
numbers. 

If an additional bound state were created by the in- 
teraction Q then the characteristic size of such a state 
should be of the order of the range of the interaction 
A (or bigger). For the interaction range, for which this 
experiment establishes a competitive limit, one obtains 
w « koX « 0.1 for A = 10 nm and w « fco'^ ~ 0.01 for 
A = 1 nm. As we mentioned previously, if such a state 
exists it would be detected in this experiment. 

Moreover, these numbers represent a lower estimation 
for the occupation numbers because they assume the 
state to be deeply bound. For a just appearing state, the 
characteristic size xq would be sufficiently bigger than 
the range of the interaction and would be close to that 
of the unperturbed states in the gravitational potential. 
Thus the population would not be almost suppressed. 



V. A LIMIT FROM EXOTIC ATOMS 

To complete our analysis in the nanometer range, let us 
note that a competitive limit follows from a recent exper- 
iment with antiprotonic atoms as well. The idea of such 
an analysis is the same as in ref. 0| . An additional force 
between a nucleus of mass M and an antiproton would 
change the spectrum of such an atom. The effective orbit 
radius rg for usually studied antiproton-nucleus atoms is 
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about a few hundreeds fm. Therefore, for the range of a 
hypothetical Yukawa force A 3> tq, an additional poten- 
tial would have a 1/r behaviour and would change the 
Coulomb potential: 



Ze^ Ze^ + mMGaa 



(31) 



The energy spectrum E„ 



2rfi 



/i of an atom with a 



reduced mass \x as well as the transition frequences 
would be modified as well: 



2mMGa, 



G 



Ze^ 



6v„ 



Thus, the constant ac is related to the shift 5vn 
transition frequencies by a formula 



Ze^ Svn 
2mMG Un 



(32) 



in the 



(33) 



An upper experimental limit on the transition frequency 
shift \5vn\ imposes an upper limit on lacj. The most 
precise measurement of the energy spectrum of antipro- 
tonic atoms was done for p'^He+ and p^He"*" atoms by the 
ASAKUSA collaboration at the Antiproton Decelerator 
at CERN J^]. In this experiment, 13 electromagnetic 
transitions between different levels of these atoms were 
investigated. No deviation was found from the values ex- 
pected within the QED calculations [liij . Equation H33|) 
allows ac to be extracted. The averaged over these 13 
transitions value of ac is equal: 



ac = (3± 10) • 10' 



27 



(34) 



and is compatible with zero, 
from this experiment is: 



A 1(7 upper limit on jad 



lacl < 1.3- 10^ 



(35) 



This value is presented by the solid horizontal line in 
Figure ^ 



VI. CONCLUSIONS AND PERSPECTIVES 

An upper limit to an additional attractive force is es- 
tablished from the measurement of quantum states of 
neutrons in the Earth's gravitational field. Relatively 
high sensitivity of the experiment 0| to a hypothetical 
additional force is due to the following factors: firstly, 
no "background" electromagnetic interactions; secondly, 
the characteristic size of the neutron wave function in the 
quantum states fits well to the range of interest for the 
short-range forces; finally, non-negligible probability to 
find neutrons (quantum-mechanical object) at distances 
much closer to the mirror than the average value of 10 
/im. 

The limit Q improves the existing contraints 9] in 
the nanometer range even if this experiment was neither 
conceived nor optimized to establish this limit. 

However, it can be easily improved in the same kind 
of experiment with some evident modifications, for in- 
stance, one can choose a mirror material (coating) with 
higher density. A more significant gain in the sensitiv- 
ity could be achieved in a dedicated neutron experiment 
which will be presented in a forthcoming paper. 
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